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Abstract. We study the Cauchy problem for the Laplace equation in a cylin- 
drical domain with data on a part of it's boundary which is a cross-section of 
the cylinder. On reducing the problem to the Cauchy problem for the wave 
equation in a complex domain and using hyperbolic theory we obtain explicit 
formulas for the solution, thus developing the classical approach of Hans Lewy 
(1927). 
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Introduction 

The question of the well-posedness of the Cauchy problem was first raised by 
Hadamard who proved in [Had23 that it is ill-posed in the case of linear second 
order elliptic equations. Hadamard's proof is based on the analytic regularity of 
linear boundary value problems. This regularity has been extended to nonlinear 
elliptic equations in [Mor58j so that Hadamard's argument also applies to general 
nonlinear elliptic equations. 

Hadamard also pointed out in jHad23j that the problem occurring in wave prop- 
agation is not at all analytic problem, but a problem with real, not necessarily 
analytic data. For general linear equations it is well known that the hyperbolic- 
ity is a necessary condition for the well-posedness of the noncharacteristic Cauchy 
problem in C°°, that is for the existence of solutions for general C°° data, cf. 
|Lax57| . |Miz61| . Moreover, for several classes of nonhyperbolic equations, explicit 
conditions on the initial data necessary for the existence of solutions were given in 
[Nis84j . For nonlinear equations, [WakOlj proves that the existence of a smooth 
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stable solution implies hyperbolicity, stability meaning that one can perturb the 
initial data and the source terms in the equations. 

The nonlinear theory yields difficult new problems, see |HF92j . jMet06| . etc. 
There are many interesting examples, for instance in multiphase fluid dynamics, 
where the equations are nor everywhere hyperbolic. As but one occurrence of 
this phenomenon, we consider Euler's equations of gas dynamics in Lagrangian 
coordinates 



mentioned in |Met06] . The system is hyperbolic, when p'{u) > 0, and elliptic, 
when p'{u) < 0. For van der Waals state laws, it happens that p is decreasing 
on an interval A mathematical example is p{u) = u{u'^ — 1). Hadamard 

argument shows that the Cauchy problem with data taking values in the elliptic 
region is ill-posed. If u{0,x) — uq{x) is real analytic near x and Uo{x) belongs to 
the elliptic interval, then any local solution is analytic, see e.g. |Mor58| . Thus, 
the initial data uq(x) must be actually analytic for the initial value problem to have 
a solution. 

It was Hans Lewy who first used hyperbolic techniques to study problems for 
elliptic equations, cf. |Lew29j . The solutions of elliptic equations with real analytic 
coefficients prove to be real analytic, and so they extend to holomorphic functions 
in a complex neighbourhood of their domain. For a holomorphic function obtained 
in this way the derivative d/dxk just amounts to the derivative d/d{iyk) where 
Zk = Xk -\- lyk are complex variables with k = 1, . . . , n. One can go to a complex 
space in only one variable, say a;„, and the change d/dxn >-> —id/dyn leads to a 
drastical modification of the characteristic variety. The Laplace equation written in 
the coordinates {x' ,Xn) with x' = {xi, . . . ,Xn-i) transforms to the wave equation 
in the coordinates {x',yn)- 

This idea is especially useful in the study of the Cauchy problem for elliptic 
equations. This problem is overdetermined even in the case of data given on an 
open part of the boundary, hence it does not admit any simple formulas for so- 
lutions, see however [Yar75| and jShl92| . Since the problem is unstable, the left 
inverse operator fails to be continuous. On the other hand, the Cauchy problem 
for hyperbolic equations is of textbook character and it admits many explicit for- 
mulas for solutions like d'Alembert, Kirchhoff, Poisson, etc. formulas, cf. |Had23) . 
Outstanding contribution to the Cauchy problem for hyperbolic equations is due 
to Leray who developed multidimensional residue theory in complex analysis to 
handle the problem, see |Ler57) . |Ler63) . etc. Having granted a solution u{x',iyn) 
of the Cauchy problem for a hyperbolic equation, how can one restore the solution 
u{x' , Xn) of the Cauchy problem for the original elliptic equation? The simple sub- 
stitution tyn 1—^ Xn does not make sense in general. For this purpose we invoke a 
formula of jCar26j which restores the values of holomorphic functions in a corner 
on the diagonal through their values on an arc connecting to faces of the corner. 
The resulting formula for the solution of an elliptic Cauchy problem includes a limit 
passage and agrees perfectly with the general observation that the character of in- 
stability in an elliptic Cauchy problem is similar to that in the problem of analytic 
continuation, cf. |Tar95) . 

As mentioned, the idea to use hyperbolic formulas for elliptic Cauchy problems 
goes back at least as far as |Lew29) . In the 1960s it was directly applied in a number 



{ 



dtu + dxv = 0, 
dxp{u) +dtv = 



(0.1) 
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of papers by Krylov, see for instance |Kry69| . In |Kry69| , an integral representation 
for hofomorphic solutions of a partial differential equation in a complex domain is 
constructed through the Cauchy data of solutions on an analytic surface. However, 
the formula does not manifest any instability of the Cauchy problem, which shows 
its local character. 

The approach we develop in this paper has the advantage of providing a large 
parameter to perturb the solution of the problem. This might give rise to a calculus 
of Cauchy problems for elliptic equations. Since these problems are unstable, no 
operator calculus similar to that including elliptic boundary values problems and 
their parametrices on compact manifolds with boundary is possible. On introducing 
a large parameter into operators we are able to describe their perturbations which 
lead to solutions. 

Let us dwell on the contents of the paper. In Section [1] we formulate the Cauchy 
problem for a second order elliptic equation in a domain X in K". The principal 
part of the equation is given by the Laplace operator while the lower order part 
may include nonlinear terms. The Cauchy data are given on a nonempty open set 
S of the boundary. Our standing assumption is that X is a, cylinder over a bounded 
domain B with smooth boundary in the space ]R"~^ of variables x' and S a smooth 
cross-section of X . 

In Section [5] we reformulate the same Cauchy problem for a hyperbolic equa- 
tion. Namely, we assume that the solution u{x' ,Xn) is a real analytic function of 
Xn G {b[x'),t{x'y) for each fixed x' e B. Then it extends to a function u{x',Zn) 
holomorphic in a narrow strip — £ <?/„<£ around the interval {b{x'),t(x')) in 
the plane of complex variable Zn — Xn + Wn- The Cauchy- Riemann equations 
force u{x' , Zn) to fulfill {d/dxn)u = —i{d/dyn)u in the strip {b{x'),t{x')) x {—e,e). 
Hence, we rewrite the original elliptic equation as a hyperbolic equation for a new 
unknown function of variables (x',?/„). Since S is the graph of some smooth func- 
tion Xn — t(x') on B, the Cauchy data transform easily for the new unknown 
function. 

In Section [3] we test our approach in the case of two variables. It is precisely the 
case treated in [Lew29| . and the approach of jLew29| does not work for n > 2. For 
rt = 2, the geometric picture is especially descriptive because the complcxification 
of X2 does not lead beyond R"^ . 

On solving the Cauchy problem for a hyperbolic equation in a conical domain in 
the space of variables {x',yn), we are left with the task of continuing the solution 
given on the base of an isosceles triangle analytically along the bisectrix of the angle 
at the vertex, for each fixed x' ^ B. To this end we invoke the classical formula of 
Carleman established precisely for this configuration, see [Car26] . Of course, the 
use of Carleman's formula is justified only for real analytic solutions of the origi- 
nal elliptic Cauchy problem. In Section |4] we give a simple proof of this formula. 
Numerical simulations with Carleman's formula failed to manifest its striking ef- 
ficiency. However, nowadays more efficient formulas of analytic continuation are 
available, cf. jAiz93| . 

In Section [5] we investigate the Cauchy problem for the inhomogeneous Laplace 
equation in the space M" of variables {x', x„) with odd n. As is shown in Section [21 
it reduces to the Cauchy problem for the inhomogeneous wave equation in the space 
of variables (a;',y„). The case n = 1 deserves a special study, for it concerns the 
initial problem for ordinary differential equations. If n = 3, the Cauchy problem 
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for the wave equation possesses a very explicit solution constructed by Poisson. For 
odd n > 5 an explicit solution formula was derived by Hadamard in |Had23| by 
his method of descent. On substituting it into Carleman's formula and changing 
integrations over y„ and a;', we get a formula for solutions of the Cauchy problem 
for harmonic functions. 

In Section |6] we restrict our attention to the Cauchy problem for the inhomoge- 
neous Laplace equation in the space R" of variables {x',Xn) with even n. By the 
above it reduces to the Cauchy problem for the inhomogeneous wave equation in 
the space of variables (x',y„). The latter Cauchy problem admits a very explicit 
solution formula due to d'Alembert in the case n — 2 and Kirchhoff in the case 
n ~ 4. For general even n the formula seems to be first published in jHad23j . We 
combine it with Carleman's formula and change the integration over ?/„ and over 
x' . This yields an explicit formula for solutions of the Cauchy problem for the 
inhomogeneous Laplace equation. To our best knowledge, this formula has never 
been published. 

In Section [7] we analyse if our approach applies to Cauchy problems for elliptic 
equations of order different from two. Yet another question under study is whether 
the method of quenching functions in the Cauchy problem for the Laplace equation 
presented in [Yar75j is actually a very particular case of formulas elaborated in this 
paper. 



1. The Cauchy problem 

Let A" be a bounded domain with piecewise smooth boundary in M". We require 
X to be of cylindrical form, i.e.. A" is a part of the cylinder B x M. intercepted 
by two surfaces yn = b{x') and 2/„ = ti^') over B, where i? is a bounded domain 
with smooth boundary in the space of variables x' = {xi, . . . , Xn^i)- For 

simphcity we assume that t{x') > b{x') for all x' G B, the case t{x') = b{x') for 
some or all x' € dB is not excluded. The Cauchy data will be posed on the top 
surface S :— {{x' ,t{x')) : x' & B} which is tacitly assumed to be real analytic, cf. 
Fig. [II 




Fig. 1. A typical domain under consideration. 
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For an elliptic second order differential operator on the closure of X the Cauchy 
data on S look like 

{u = uo on 5, 
du 
— = ui on S, 
ov 

where v is the outward unit normal vector at S. Obviously, v — Vgi/|Vp| where 
g = Xn - t{x'). 

Lemma 1.1. If u is a smooth function near S satisfying u — uq on S, then 

9^ a/|V:,4|2 + 1 V dXn I 

on S. 

Proof. This is an easy exercise. □ 

Consider a nonlinear second order partial differential equation Au — f{x, u, Vu) 
in X, where f{x, u,p) is a real analytic function on X xRx R". By Lemma [TTTl the 
Cauchy problem for solutions of this equation with data on S can be formulated in 
the following way. Given functions uq and ui on S, find a function u in X smooth 
up to S which satisfies 

Au — /(a;, li, Vu) in X, 
u — Uf) on S, (1-1) 

u'^^^ = ui on S. 

Lemma 1.2. There is at most one real analytic function u in X U S which is a 
solution of 

Proof. Let ui and U2 be two real analytic functions in A" U 5 satisfying (jl.ip . Set 
u = ui — U2, then u is real analytic uy X \J S and vanishes up to the order 2 on 5. 
Hence it follows that Au = f{x, ui, Vwi) — /(x, U2, VU2) vanishes on S. Since A is 
a second order elliptic operator, we readily deduce that u'^^^^ = on 5, and so u 
vanishes up to order 3 on iS. Hence it follows that Au vanishes up to order 2 on 5, 
and so (d/dxn)^u = on iS. Arguing in this way, we conclude that u vanishes up 
to the infinite order on S. Since u is real analytic in A" U 5, we get u = in A", as 
desired. □ 



2. Hyperbolic reduction 

Assume that u is a real analytic function in XUS which satisfies . Then, for 
each fixed x' G B, the function u{x' , Xn) can be extended to a holomorphic function 
u{x' , Xn+iyn) in some complex neighbourhood of the interval {b{x'), t{x')]. Without 
loss of generality we can assume that this neighbourhood is a triangle T{x') in the 
complex plane z„ = x„ + iy„ with vertexes at b{x') and t{x') =p le, where e > 
depends on x' . We write U{x' ,Xn,yn) for the extended function, so that u(x) just 
amounts to J7(x', a;„, 0). 

Since u(x',Zn) is holomorphic in a complex neighbourhood of {b{x'),t{x')], it 
follows from the Cauchy-Riemann equations that 

(sx") " ( ^ U{x',Xn,yn) 
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for all j — 1,2,.. .. Therefore, the Cauchy problem for u transforms to the 

problem 

r A,,u-u;;^y^^ = fix',zn,u,v,,u,-iu;^j, if x' e s, z„ e t(x'), 

< U{x',Xn,0) = Uo{x',Zn), if x' E B , Zn ^ t{x') , 

[ Uy {x',Xn,0) = lUi{x',Z„), if x' e B, Zn ^ t{x'), 

(2.1) 

relative to the new unknown function C/(a;', x„, 2/„). 

Hardly can p.ip be specified within Cauchy problems for second order differential 
equations, for the number of independent variables is n + 1 while the Cauchy data 
are given on a surface of dimension n — 1. Since the differential equation in (12. ip 
does not contain the derivative U'j. , it is easy to deduce that the smooth solution to 
this problem is by no means unique. This no longer holds true for the holomorphic 
solution because of uniqueness theorems for holomorphic functions. Moreover, if 
U{x' , Xn, Un) is holomorphic in Zn = Xn + Wn, then the differential equation in ()2.1|) 
is satisfied for all x' G B and z„ S T{x') provided it is fulfilled for all x' E B and 
Zn = t{x') + tyn with |y„| < e. 

Thus, when one looks for a holomorphic solution to (|2.ip . this problem actually 
reduces to the Cauchy problem for a quasilinear hyperbolic equation in the space 
of variables (a;',y„), whose principal part is given by the wave operator. More 
precisely, 

U^^y^ = A,,U - f{x',Xn+iVn,U,V^,U,^iU'yJ, ff x'^B, 

\yn\ < e(a;'), 

U{x',x„,0) = UQ{x',Xn), if x' e B, 

Uy {x',x„,0) = iui{x',Xn), if x' e B, 

(2.2) 

where the variable Xn is thought of as a parameter which runs over the interval 
{b{x'), t{x')). We are actually interested in the solution of this problem correspond- 
ing to the special choice Xn = t{x') of the parameter. In other words, we study 
problem (|2.2p on the hypersurface a;„ = t{x') in the space of variables {x, yn), the 
Cauchy data being given on the intersection of the hypersurface with the hyperplane 

{yn = 0}. 

When passing to the Cauchy problem on the hypersurface x„ = t{x') in R""*"^, 
one should interpret equations (|2.2p adequately in accordance with the presence of 
parameter x„. Namely, each equations has to be fulfilled together with all deriva- 
tives in Xn on x„ — t{x'). 

Lemma 2.1. There is at most one function [/(x', x„, y„) in a neighbourhood ofS, 
which is real analytic in yn at yn = and satisfies h2. 2\) with Xn — t{x'). 

Proof. Let Ui and U2 be two functions in a neighbourhood of S, which are real 
analytic in yn at yn = and satisfy (|2.2|) with x„ = t{x'). In the coordinates 
{x',Xn,yn) the surface S is given as intersection of two hypersurfaces x„ — t{x'), 
where x' G B, and ?/„ = 0. Set U — Ui~U2, then U is real analytic in y„ at j/„ = 0. 
We shall have established the lemma if we prove that each derivative (9/9j/„)-' C/ 
with j = 0, 1, . . . vanishes for a;„ = and ?/„ — 0. For j — 0,1 this follows 

immediately from the conditions which Ui and U2 fulfil on S. For j < 2 this follows 
from the differential equation in p.2|) by induction. We check it only for the initial 
value j = 2, for the induction step is verified in much the same way. From (|2.2|) we 
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get 

provided that a;„ — t(x'). 

Since (9/9a;„)^ (C/i — U2) = for a;„ = t{x'), y„ = 0, and all j = 0,1,..., it 
follows that 

Ul,^{x',t{x'),{)) = (;7i(x',t(x'),0)4-C/i,,„(a;',t(x'),0)t;^(a;') 

= {U2{x' ,t{^'),{))t,-U2,.S^' ,t{x'),Q)t',^{x') 

- C/L.(2:',i(a;'),0) 

for each /c = l,...,n— 1. Moreover, we get 

5^;[/i=9,";C/2 (2.3) 

on the surface a;„ = t{x'), yn = for all multi-indices a' — (ai, . . . ,a„_i). This 
yields readily Ax'Ui — Ax'U2 for x„ — t{x') and ?/„ = 0. Substituting these 
equalities into the formula for Uy^y^ we obtain Uy^y^{x' , t{x'), 0) = for all x' G B, 
as desired. □ 

Note that equalities (|2.3p generahse to d"dy^^^Ui = d^dy^^^U2 for x„ = t{x'), 
yn = 0, and all multi-indices a = (ai, . . . ,a„) and a„+i = 0, 1, . . ., as is easy to 
check. 

We have thus reduced the Cauchy problem for the Laplace equation perturbed 
by nonlinear terms of order < 1 to the Cauchy problem for the wave equation 
perturbed in the same way. The reduction is justified as long as the solution under 
study is real analytic in x„. 

Perhaps the reduction does not make sense in the case n = 1, for it leads to no 
simplification. 

3. The planar case 

To test the hyperbolic reduction of Section [2l we consider the case n = 2 in 
detail, assuming / to depend on x £ X LI S only. 

Let A' be a strip domain in consisting of all x — (xi, X2), such that xi E (a, b) 
and b{xi) < X2 < t{xi), where (a, 6) is a bounded interval in R and b, t are 
smooth functions of xi E (a, &). Write B {a,b) and denote by S the curve 
{{xi,t{xi)) : xi e (a, 6)} which is a part of dX. We focus on the Cauchy problem 
for the inhomogeneous Laplace equation given by (|l.ip . When looking for a solution 
u of this problem which extends to a holomorphic function u(xi , 22) of Z2 = a;2 -f iy2 
in a neighbourhood of {(2:2,0) : X2 € (b(xi),t(xi)]}, for each fixed xi £ {a,b), we 
arrive at 

Uy^y2 = U^.xi- f{xi,X2 + iy2), if xiE{a,b), 

\y2\<£(xi), 

U{xi,X2,0) = Uq{xi,X2), if Xi E ia,b), 

Uy^{xi,X2,0) = tUi{xi,X2), if XiE(a,b), 

which is a Cauchy problem for the inhomogeneous wave equation with parameter X2 
relative to the unknown function U{xi,X2,y2) = u{xi,X2 + cf. (|2.2p . We are 
actually interested in finding a function U which satisfies (|3.ip only on the surface 
X2 = t{xi), see Fig.m 
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Fig. 2. The case n = 2. 



It is an easy exercise to verify that the function 

(G/)(xi,a;2,2/2) = --r / dy'^ f{x[,X2+i(y2-y'2))dx[ 

^ Jo Jxi-y'^ 

satisfies the inhomogeneous wave equation and homogeneous (i.e., corresponding to 
1*0 = ^1= 0) initial conditions in (jS.ip . On the hand, d'Alembert's formula gives a 
function satisfying the homogeneous (i.e., corresponding to / = 0) wave equation 
and the inhomogeneous initial conditions in p.ip . see |CH68[ Ch. I, § 7.1]. In fact, 
this is 

^/ X uo{xi+y2,X2) +uo{xi-y2,X2) I ri+y^ 

P{uo,Ui){xi,X2,y2) ^ ^ 1- ;r / Ui(Xi,X2)dXi, 

(3.2) 

where the right-hand side is well defined for all (xi, X2, 2/2) satisfying Xi+y2 G (a, b) 
and Xi — y2 G (a, b). The pairs {xi , 2/2) with this property form two cones in the 
plane, being the set of all {xi,y2), such that xi S {a,b) and ±2/2 € [0,e(a;i)), 
where 



e(a;i) 



Xl 



2 

Thus, given any twice differentiable function uq(xi,X2), differentiable function 
ui(xi,X2) of Xl € (a, 6) and any differentiable function f{xi,Z2) of both variables, 
the formula 

U = Gf + Piuo,ui) 

yields a solution to the Cauchy problem (|3.1I) for all values of parameter X2 that 
do not lead beyond the domains of uq, ui and /. Had we known uo(xi,X2) and 
ui{xi,X2) for all values X2 G {b(xi),t{xi)], then the first initial condition of (13. ip 
would give U{xi , X2 , 0) — uq {xi , 2:2) and so the solution to the Cauchy problem (II. ip 
by it(x) = uo{xi, 2:2). This just recovers the reduction but is not of use to solve the 
original Cauchy problem. However, on substituting X2 = t{xi) into U{xi,X2,y2) 
we obtain 

I ry2 rxi+y'2 

u{xi,t{xi) +iy2) ^ I dy'2 f{x[,t{xi)+t{y2-y'2))dx[ 



2 



UQ{xi+y2,t{xi)) +uo{xi-y2,t{xi)) i f ^^^^ , , 
+ -/ ui{xi,t{xi))dx^ 

(3.3) 
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for all xi E {a,b) and \y2\ < e{xi). Note that {x'i,t{xi)) fails to lie on the curve S 
for all x[ S [xi — 1/2, xi + 2/2] unless t{xi) is constant. Therefore, u{xi,t{xi) + iy2) 
is determined by the Cauchy data of u in some neighbourhood of S. This forces us 
once again to confine ourselves with solutions which are real analytic in the variable 

X2. 

For fixed Xi € (a, &), formula p.3|) gives the restriction of the function u{xi, Z2), 
holomorphic in Z2 in the triangle with vertexes at b{xi) and t{xi) =F is{xi), to the 
side t{xi) + t[—e{xi),e{xi)] of the triangle. This limits application of hyperbolic 
theory. Our next objective is to continue the function from the side of the triangle 
analytically along the bisectrix of the angle at b{xi). This is a problem of analytic 
continuation. 



4. Carleman formula 



Let be a domain in the complex plane C of variable z bounded by lines BO and 
OA and by a smooth curve c = AB lying inside the angle BOA. Write ZBOA — an 
with < a < 2. 

Choose the univalent branch of the analytic function in the complex plane 
with a slit along the ray aigw = tt, which takes the value 1 at it; = 1. 

Lemma 4.1. If u is a holomorphic function in D continuous up to the boundary, 
then 



u(z) = lim - — / u(0 exp A'^ 

N->-oo ZTTl 



^ - Co 



c- 



holds for any point z E D on the bisectrix of the angle BOA, where Co is o, complex 
number corresponding to the vertex O of the angle. 

This formula is due to Carleman |Car26| . To our best knowledge it was the 
first formula of analytic continuation using the idea of quenching function. Since 
that time such formulas in complex analysis and elliptic theory are called Carleman 
formulas, see |Aiz93| . |Tar95) . 



Proof. Fix any z £ D lying on the bisectrix of the angle BOA. For N 
we apply the Cauchy integral formula to the function 

'C-Co\i/" 



1,2,.. 



i(C) exp7v((- 



1 



- z - Co ' 

which is holomorphic in D and continuous in the closure of D. Since its value at 



C = z is u{z), we get 
u{z) = 



— / u(C) expA^ 



1 

27ri 
1 

271^ 



- 1 



dD\c 



C-Co\i/" 



dc 



-Co 



(4.1) 



IfC edD\c, then 



Vz-Co/ 





c- 


Co 




z — 


Co 



1/c 



exp 



± 



C-Co 

z-Co 



1/c 
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i^-^ — - Ij equals e^^. Letting N 
(|4.ip establishes the lemma. □ 

Having disposed of this preliminary step, we now turn to the problem of analytic 
continuation we have encountered in Section [3] We apply Lemma 14.11 in the plane 
of complex variable Z2 — X2 + iy2- Given any fixed Xi G (a, &), we take the triangle 
T{xi) with vertexes O :— h{xi) and A :— t{xi) — i£{xi), B := t{xi) + ie{xi) as D, 
cf. Fig. [21 




Fig. 3. Recovering a holomorphic function. 



a = -arctanf^4i^ 



In this case 

2 

TT \t{xi) — b{xi)' 

depends on xi and the bisectrix of the angle BOA coincides with the real axis. The 
solution u{xi,Z2) is given on the edge AB and we are aimed at reconstructing it in 
the interval (b{xi),t{xi)). 

Theorem 4.2. Let n = 2. For each solution u of the Cauchy problem in X 

which is real analytic up to S, the formula 



u{x) ~ lim 



^ I TT, .( \ \ Ari^i^ ^(^i)-K^i)+»y2 ^ g A dy2 

— / [/(a;i,t(a;i ,y2 expiV — -1 TTT 

'21: J W X2 ~ b[Xi) ) ) t(x{) 



N^oo2lT J W X2 ~ b[Xi) / t{xi) — X2+iy2 

-e(xi) 

holds for all X €z X . 

Proof. This follows immediately from Lemma 14.11 and formula (jS.Sp giving an ex- 
plicit continuation of the solution u{xi , X2) along S to the plane of complex variable 

Z2 = X2 + iy2 ■ □ 

This formula is especially simple if 5 is a segment X2 ~ to, i.e. the graph of a 
constant function t{xi) = to of xi € {a,b). If moreover / = then formula p.3p 
transforms to 

jrr , ^ uo{xi+y2,to) +uo{xi-y2,to) i f'^^+y^ 1 1 + \^ 1 
U{xi,ta,y2) = 5 h-/ ui{x^,ta)dx^ 
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for all xi G ia,b) and |j/2| < £{xi). Substituting this into the formula of Theorem 
4.21 we get 

u{x) = lim / u{x\,tf))'R Kn{xitX2,Xi — x'l) dx'i 



N- 



lim 



(a;i,io)( / KNixi,X2,y2)dy2^dx[, 



N^oo J dx2 

xi-e{xx) \x'-^^-xi\ 

(4.2) 



where 



t{x') ~ b{x') + 3 _ ^ 



, / \ ^ ^ ^ Xn b[x'} 

Kn{x ,Xn,yr ^ - 



2tT t{x') - Xn + lyn 

Formula (|4.2p can be regarded as an elliptic analogue of the d'Alembert formula 
for the wave equation. 

Note that nowadays there are many explicit formulas of analytic continuation 
which are simpler than the original formula of |Car26| . We refer the reader to 
|Aiz93) . 

5. POISSON FORMULA 

In this section we discuss the case n = 3 in detail, assuming the function / to 
depend on x G A" U 5 only. The Cauchy problems for the inhomogeneous Laplace 
equation reduces to the Cauchy problem for the inhomogeneous wave equation. 
This latter reads 

Kv. = A^,-f(x\x^ + iy:i), if x'eB, 

\yz\<e{x'), , . 

Uix',X3,0) = uo{x',X3), if x'eB, ' 

^Uy^{x', X3,0) = iui{x',X3), if x'eB, 

X3 being thought of as parameter. We are aimed at finding a function U which 
fulfills (|5.ip on the surface 2:3 = t(x'). 

The advantage of the reduction lies in the fact that the Cauchy problem for 
hyperbolic equations is well posed in the class of smooth functions. For n — 3, 
there is an explicit formula for its solution due to Poisson, see |CH68| Ch. Ill, 
§ 6.5]. More precisely, 

jr, I 1 /j / / f{x",X3 + i{y3-y'^)) „ 
U(x ,X3,y3) ^ -— dy^ / ======- dx 



c"~x'\<\y'^\ 

d sgnj/s f uo{x",X3) ^ sgny3 f mi{x",X3) 



\x" -x'\<\y3\ \x" -x'\<\y3\ 

(5.2) 

for all x' £ B and jya] < e{x'). 

For formula (15. 2p to make sense it is certainly required that, for any 2/3, the ball 
\x" — x'\ < 1 2/3 1 would belong to the domain B in M"r^, where the Cauchy data 



12 



D. FEDCHENKO AND N. TARKHANOV 



uo{x\xn) and ui{x',Xn) are given. Since ys varies in the interval {—e{x'),e{x')), 
we get readily the formula e{x') — d{x' , dB), the distance from x' to the boundary 

of B, cf. Fig. m 



d{x',dB) 




Fig. 4. Reduction to imaginary cones. 



Theorem 5.1. Let n — 3. For each solution u of the Cauchy problem in X 

which is real analytic up to S , the formula 

e(x') 

2 / e(x') 



-1 



JV^oo 27r ./ \ ' \ JTa^Jj t- \^\^ _ 5(^2;') J J l(^x') — X3 + iyj, 



holds for all X £ X , where a = — arctan ( — — - 

TT \tx') 



t{x')-b{x')- 

Proof. This is a direct consequence of Lemma 14.11 and formula (|5.2p which gives 
an explicit continuation of the solution u(x',X3) along S to the plane of complex 
variable — x^ + ly^. □ 

On substituting (|5.2p into the Carleman formula of Theorem 15.11 we arrive at 
an explicit formula for solutions of the Cauchy problem for the inhomogeneous 
Laplace equation. The computations are cumbersome, and so we confine ourselves 
with the case / = 0, as in (|4.2I) . By the very construction of the Carleman kernel, 
Kn{x' ,x^,e{x')) tends to zero as — >■ oo, for any x' £ B and 0:3 £ {b{x'),t{x'j). 
Hence 

d_ 

dys^ 



u{x) 



lim 

AT— >oo 



u{x",t{x'j) 



^^?fiKN{x',X3,y3) 



\x" —x'\<e{x') 



— lim 

A->oo 



\x" —x' I 

e(x') 



\x" —x'\<e{x'') 



\x" —x' I 



^ ^yi-\x" -x'\^ 
1 '^Kf^{x' ,x-i,y-i) 



dyz dx' 



dys ) dx" 



(5.3) 



for all X £ X. 

Formula (|5.3p can be thought of as an elliptic analogue of the Poisson formula 
for the wave equation. 
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6. KiRCHHOFF FORMULA 

The solution of the Cauchy problem for the wave equation bears certain structure 
which changes in odd and even dimensions. For this reason we consider also the case 
n = 4 in detail. The corresponding formula for solutions of the Cauchy problem 
for the wave equations is known as the Kirchhoff formula, see |CH68| Ch. Ill, § 6.4] 
and elsewhere. 

By the above, the Cauchy problem for the Laplace equation in a cylindrical 
domain X (Z&f^ reduced to 

U^\y^ = A,^, - /(x',X4 + lyi), if x' e B, 

\yi\<e{x'), 

f/(a;',X4,0) = Uo{x',X4,), if x' e B, ' 

^Uy^i^', X4,0) = iUi{x',X4), if x'eB, 

where x' = {xi,X2,X3) varies in a domain B C M^, £{x') stands for the distance 
from x' € B to the boundary of i3, and xs is thought of as parameter in (b{x'), t{x')]. 
The Cauchy data uq and ui are in C^{B) and C^{B), respectively. The Kirchhoff 
formula gives 



An _ 

\x" —x' \ <\y4\ 

+ T{ — -; / uo(x" ,Xi)dcr(x") + / tui(x" , X4,)da(x") 

dyi 47ry4 J Airy^ J 

\x" -x'\ = \y4\ \x"-x'\ = \y4\ 

(6.2) 

for all x' E B and \y4\ < e{x'). 

The substitution Xi = t{x') into U gives the restriction of the function U , holo- 
niorphic in Z4 = x^+iy^^ to the edge £(2:'), £(a;')] of the triangle r(x') C C, 

where U is holomorphic. Using Carleman's formula of Lemma 14.11 we arrive at a 
formula for u{x) similar to that of Theorem l5.ll It reads in much the same way, with 
x^ and j/3 replaced by X4 and 1/4, respectively. For short we restrict our attention 
to a formula like (|57 



Corollary 6.1. Let n = 4. For each solution u of the Cauchy problem with 
f = Q in X , which is real analytic up to S , we get 

r 1 (-S-^Kn^{x\x4,\x" - x'\) 

u{x) = ~ lim / u{x",t{x'j)—-^ rv, 7, dx" 

^ ' N^oc J ^ ' ^ " 2tt \x" - x'\ 

\x" —x'\<e{x') 

_ lim f ^"fj-" t(x']] ^ ^^KNix',X4Ax" -x'\) 

J dx4^ ' ^ ''27T \X"~X'\ 

\x" —X' I <6{x') 

(6.3) 

for all X E X . 



Proof. The proof is quite elementary although cumbersome. We first substitute the 
integral of mq on the left-hand side of (|6.2p into Carleman's formula. Integration 
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by parts yields 

e{x') 



J '^iiny' J uoix",t{x'))da{x")jKN(x',X4,y4)dy4, 

\x"-x'\ = \y4,\ 

i.ATTy' J uo{x" ,t{x'))da{x")^KN{x' ,Xi,yi 

i-x'\ = \yi\ 

1 f „,\ d 



x"-x'\ = \yi\ 

e(x') 

, uo{x",t{x'))da{x"))——KNix',X4,y4)dy4:. 
Airyi J / dyi 

-e{x') \x" -x'\ = \y4\ 

The first integral on the right-hand side is equal to 

(2^^ / uoix",t{x'))daix")^^KNix',Xi,e{x')), 

\x" — x' |— e(a:') 

which vanishes as iV — !■ 00 by the construction of the kernel ifAr(a;', X4, e(a;')). 
Indeed, the point t(x') + leix') belongs to the top leg of the angle BOA, and X4 to 
its bisectrix. 

Furthermore, we write the second integral on the right-hand side as the sum of 
two integrals. The first integral is over y^ e (— e(a;'),0) and the second one over 
2/4 G (0,e(a;')). In the second integral we change the variable by 2/4 !->■ —2/4, and 
then evaluate the sum, obtaining 

e{x') 

/ uq{x" ,t{x'))da{x"))- — Kn{x' ,X4:,yi) dyi 

. 47rt/4 J J dy4 

-e(a:') \x" -x'\ = \y4\ 

^^^'^ 1 d 

(- / uo{x",t{x'))da{x"))——diKN{x\x4,y4)dy4. 

\ 27r2/4 J ) dy4 

\x" ~x'\ = \yi\ 

Since dx" = da{x")dy4, we deduce from Fubini's theorem that the latter integral 
just amounts to 

(^}fjr„)(i'.i4,k"-i'l) 

\x" — x' \ <e(x') 

as desired. 

The same (even easier) reasoning applies when one substitutes the integral of ui 
on the left-hand side of (|6.2p into Carleman's formula. The details are left to the 
reader. □ 



Formula (j6.3l) is an exposition of Kirchhoff 's formula for the wave equation in 
the context of elliptic theory. We have already mentioned another interpretation of 
Kirchhoff's formula in |Kry69l . Unfortunately, we could not understand this latter 
paper. 
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7. Concluding remarks 

The developed method of analytic continuation in the plane of complex variable 
Zn = Xn + Jj/n Still works if the Cauchy problem under study is nonlinear. Having 
granted a holomorphic solution f/(x', x„, j/„) to the Cauchy problem (|2.2p on the 
surface a;„ = t{x'), we use Carleman's formula to extend U to all of X. The 
extension looks like 

u{x) ^ lim / U{x\t{x'),yn) KNix\xn,yn)dyn (7.1) 

for all x G X. 

Formula (|7.1|) allows one to construct explicit formulas similar to (|4.2p . (|5.3|) 
and (|6.3|) for arbitrary n. To this end one uses classical formulas for the solution of 
the Cauchy problem for a second order hyperbolic equation by the descent method 
of Hadamard, cf. jHad23j . |CH68[ Ch. VI, §. 5.2]. We were rather interested in 
equations of mathematical physics. 

The simplest formula is obtained for even n > 4, thus generalising Kirchhoff's 
formula If uo £ C("+2)/2(5) and ui e C"/2(5)^ then every solution u of ^ 

with / = represents by 

u(x) — lim / dx" 



\x" — x'\<e{x') 



( n*-i9 ((^~ — ~) ' yn'SiKN){x',x„,\x"~x'\) 



fT„_il-3 •. . .• (ri-3) 



o, , iN5f-i9 I it; ' Q?ifjv (a;',x„,|a;"-x'|) 

' cr„_il-3 •. . .• (n— 3) \x" — x'\ 

(7.2) 

for all X £ X, where cr„^i stands for the area of the (n — 2) -dimensional unit sphere 
in R"~^. We used here an exotic designation for the integral by purely technical 
reasons. 



Remark 7.1. Formula ()7.2p has much in common with the familiar formula of 
|Yar75j . 

The method of proof carries over to right-hand sides f{x, u, Vu) which are afRne 
functions of u and Vu. This is the case, e.g., for the Helmholtz equation, cf. jCII68[ 
Ch. VI, §. 5.7]. 

Another class of equations which may be handled in much the same way consists 
of those of the form 

^" + <„a;„ = fix), 

where A is a linear differential operator containing at most the derivative u'^ but 
no higher order derivatives in a;„, see [CH68| Ch. Ill, § 6.4]. 
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